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ABSTRACT 

The  investigation  dealt  with  the  numerical  and  experimental 
investigation  of  the  effect  of  porosity  on  the  vortex  shedding  from  a 
cambered  plate.  The  effect  of  porosity  was  incorporated  into  the  anal- 
ysis through  use  of  the  experimentally  obtained  pressure-drop  data. 
The  results  have  shown  that  the  porosity  of  the  existing  materials  used 
in  parachute  canopies  is  not  sufficient  to  eliminate  the  negative  pres- 
sure gradient.  Canopies  with  nonuniform  porosities  will  have  to  be 
designed  to  overcome  the  collapse  phenomenon  resulting  from  the 
unfavorable  pressure  gradients. 
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I.    INTRODUCTION 

The  current  parachute-canopy  design  methodology  is  mainly 
empirical.  It  is  based  on  either  model  tests  in  wind  tunnels  (see,  e.g., 
Klimas  and  Rogers,  1977;  Pepper,  1986;  Rychnovsky,  1977)  or  on  the 
deployment  of  full-scale  parachutes  with  test  loads.  However,  under 
certain  conditions,  parachutes  have  been  observed  to  collapse  during 
periods  of  rapid  deceleration,  subsequent  to  the  initial  inflation  of  the 
canopy  with  a  resulting  loss  of  drag  and  ultimately  the  payload. 

Studies  began  with  a  two-dimensional  model  of  flow  about  a  rigid, 
nonporous  camber.  Even  though  highly  idealized,  the  purpose  of  the 
investigation  was  not  a  solution  of  an  immediate  practical  problem  but 
to  understand  the  physics  of  the  problem.  With  that  in  mind,  the  dis- 
crete vortex  model  (DVM)  (see,  e.g.,  Chorin,  1973;  Sarpkaya,  1975) 
was  used  to  model  the  flow  past  bodies  with  included  angles  of  120, 
180,  and  240  degrees,  and  was  analyzed  by  Mostafa  (1987)  and  Munz 
(1987).  These  investigations  have  shown  that  the  drag  force  is  pre- 
dicted fairly  well  during  the  steady-state  period  and  up  to  the  middle 
of  the  deceleration  period.  Beyond  that  time,  the  predicted  forces 
tended  to  be  relatively  larger  than  those  measured  experimentally. 

Nevertheless,  the  forces  measured  during  the  final  stages  of 
deceleration  were  very  small,  so  the  degree  of  exactness  was  not  a 
concern.  Moreover,  the  reasons  leading  to  the  collapse  (backwash  of 
the  vortices)  occurred  shortly  after  the  start  of  deceleration,  where 


the  predicted  and  measured  forces  compared  well.  Previous 
investigations  at  the  Naval  Postgraduate  School  have  concluded  that 
the  development  of  negative  differential  pressures  in  the  central 
region  of  the  camber  is  the  likely  cause  of  parachute  collapse  at  high 
deceleration  rates  (see,  e.g.,  Munz,  1987). 

In  view  of  the  foregoing,  it  became  clear  that  the  flow  through  the 
porous  camber,  if  managed  properly,  can  help  to  minimize  or  elimi- 
nate the  collapse  phenomenon  by  eliminating  the  negative  differential 
pressure  across  certain  sections  of  the  camber.  The  present 
investigation  was  undertaken  to  study  in  detail  the  effect  of  porosity  on 
a  two-dimensional  rigid  camber  in  order  to  bring  the  theoretical 
problem  a  step  closer  to  reality.  In  the  following,  the  fundamental 
concepts  used  in  the  introduction  of  flow  through  a  porous  wall  and 
the  results  obtained  with  various  methods  and  porosities  are  described 
in  detail. 


H.    EFFECTS  OF  POROSITY 

A.      REVIEW  OF  THE  FLOW  THROUGH  POROUS  SCREENS 

A  parachute  canopy  is  a  porous  material.  The  amount  of  flow  that 
goes  through  such  a  material  depends  on  the  percentage  area  opening, 
the  shape  and  size  of  the  threads,  the  moisture  content  of  the  mate- 
rial, the  differential  pressure  across  the  canopy,  the  angle  of  approach 
of  the  ambient  flow,  the  magnitude  and  direction  of  the  flow  on  the 
downstream  side  of  the  canopy,  and  the  degree  of  time-dependence  of 
the  upstream  and  downstream  flows  (i.e.,  the  vortex  motion  down- 
stream of  a  parachute). 

In  order  to  incorporate  the  effect  of  porosity  into  the  analysis,  one 
had  to  know  at  least  the  flow  rate  through  the  canopy  as  a  function  of 
differential  pressure,  assuming  that  all  of  the  other  parameters  cited 
above  were  of  secondary  importance,  or  that  their  effects  may  be 
incorporated  into  a  flow  coefficient  relating  differential  pressure  to 
the  normal  velocity  of  the  flow  rate.  This  seemingly  simple  search  for 
a  reliable  (theoretical  and /or  experimental)  equation  turned  out  to  be 
anything  but  simple.  The  following  is  a  brief  description  of  the  previ- 
ous investigations  of  the  flow  past  screens. 

A  fairly  detailed  review  of  the  flow  through  screens  is  given  by 
Laws  and  Livesy  (1978).  There  are  a  number  of  investigations  (see, 
e.g.,  Baines  and  Peterson,  1951;  Turner,  1969;  Reynolds,  1969; 
Graham,  1976)  which  dealt  with  the  flow  past  fairly  porous  screens  for 


the  purpose  of  understanding,  managing,  or  generating  isotropic  tur- 
bulence. Often  the  resistance  through  the  screen  was  of  minor 
importance.  Others  (see,  e.g.,  Taylor,  1944;  Taylor  and  Davies,  1944; 
Wieghardt,  1953)  who  investigated  the  flow  through  infinite  screens 
expressed  the  resistance  as 

Ap  =  ipKU2 

Koo  and  James  (1973)  and  Cumberbatch  (1982)  investigated  the  flow 
past  finite  screens  set  in  infinite  flows,  but  the  effects  of  body  shape, 
separation,  vortex  formation,  flow  magnitude  and  direction  down- 
stream of  the  body,  and  deceleration  of  the  flow  were  not  studied.  It 
became  clear  that  there  was  an  enormous  amount  of  research  done  on 
the  flow  through  screens  and  yet  there  was  very  little  information 
which  could  be  extracted  from  them  for  use  in  this  study  of  a  deceler- 
ating parachute.  Payne  (1978)  reviewed  the  data  from  many 
researchers  in  an  effort  to  relate  the  geometric  porosity  (open  area 
ratio)  to  the  differential  pressure,  and  developed  an  equation  in  terms 
of  dynamic  pressure  loss  term  and  a  viscous  loss  term, 

Ap  =  K1U2-f  K2U 

However,  this  study  involved  many  different  cloth  types,  not  just 
parachute  cloth,  and  included  no  evaluations  of  the  effects  of  flow 
direction.  This  prompted  the  undertaking  of  an  experimental  investi- 
gation to  establish  the  proper  relationship  between  the  pressure  drop, 


the  flow  rate,  and  the  flow  direction  on  the  actual  samples  of  canopy 
materials  used  in  parachutes. 

B.      PARACHUTE  CLOTH  EXPERIMENTS 

1.     Experimental  Apparatus  and  Procedure 

The  purpose  of  the  experiments  was  twofold:  to  find  the  dif- 
ferential pressure  versus  flow  rate  relationship,  and  to  find  the 
change,  if  any,  of  this  relationship  with  the  downstream  flow  direction. 
For  this  purpose,  an  experimental  apparatus  was  conceived,  designed, 
and  manufactured  as  shown  in  Figures  1  and  2.  The  flow  rates  at 
points  A  and  B  were  measured  with  the  standard  flow  meters,  and  the 
static  pressures  upstream  of  the  screen  were  recorded  for  each  flow 
rate  through  the  screen.  Static  pressure  downstream  of  the  screen 
was  assumed  to  be  atmospheric. 

Observations  have  shown  that  the  canopy  sample  did  not 
bulge  measurably  during  the  experiments,  even  at  the  differential 
pressures  encountered  by  a  parachute  beginning  to  descend  at  a  speed 
of  about  600  miles  per  hour.  This  was  primarily  due  to  the  proper 
construction  and  sandwiching  of  the  circumference  of  the  canopy 
sample.  However,  this  is  not  meant  to  exclude  the  possibility  that  at 
high  pressures,  the  slight  cambering  of  the  sample  may  have  stretched 
it  and  thereby  increased  its  porosity  slightly,  leading  to  a  smaller  dif- 
ferential pressure  for  a  given  flow  rate  than  that  had  it  not  been 
stretched.  In  fact,  this  was  checked  by  conducting  the  experiments 
first  with  increasing  flow  rate  and  then  decreasing  the  flow  rate.    No 


Figure  1.   Parachute  Cloth  Differential  Pressure 
Measurement  Apparatus,  Front  View 


Figure  2.   Parachute  Cloth  Differential  Pressure 
Measurement  Apparatus,  Three-Quarters  View 


attempt  was  made  to  wait  a  long  time  after  the  completion  of  the 
experiment  with  ascending  pressures  to  see  if  the  material  would 
return  to  its  original  unstretched  shape,  or  whether  fibers  used  in 
construction  of  the  canopies  develop  a  permanent  set  beyond  certain 
differential  pressures  and  retain  their  deformed  shape.  This  leads  to 
the  possibility  that  a  once-used  parachute  may  fall  faster  than  a  never- 
used  one. 

Of  the  seven  available  canopy  materials,  four  samples  were 
chosen  and  tested  to  provide  a  range  of  data.  The  characteristics  of 
the  four  samples  are  shown  in  Table  1.  The  identifying  letters  are 
used  for  later  reference. 

TABLE  1 
CHARACTERISTICS  OF  PARACHUTE  CLOTH  SAMPLES 


Sample 
Letter 

Material 

Weight 

Military 
Specifications 

Number 

A 

Nylon 

1.1  oz. 

— 

B 

Nylon 

4.75  oz. 

MIL-C-8021D 
Type  I 

C 

Nylon 

7.0  oz. 

MIL-C-8021D 
Type  II 

D 

Kevlar 

3.0  oz. 

— 

During  the  course  of  the  experiments,  the  flow  on  the  down- 
stream face  of  the  sample  was  systematically  changed  to  see  if  the 
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differential  pressure  across  the  sample  would  increase  for  a  given  flow 
rate  through  the  sample. 
2.     Results 

The  results  of  tests  on  the  effect  of  flow  on  the  downstream 
face  of  the  samples  have  emphatically  and  somewhat  surprisingly 
shown  that  the  flow  through  the  screens  investigated  is  not  affected  by 
the  presence  or  absence  of  the  downstream  flow.  This  apparently 
simple  finding  is  of  extreme  mathematical  and  practical  importance. 
It  simplifies  the  analysis  and  the  possible  management  of  porosity  to 
eliminate  parachute  collapse. 

The  results  obtained  with  the  four  representative  canopy 
samples  are  shown  in  Figures  3  and  4. 

Figure  3  shows  the  variation  of  differential  pressure  as  a  func- 

1  2 

tion  of  the  dynamic  head,  j pVN  ,  where  Vn  is  the  velocity  through  the 
canopy  (j)  .  It  is  clear  that  the  value  of  differential  pressure  varies 
with  the  weight  of  the  sample  and  manufacturing.  For  samples  B  and 
C,  the  military  specifications  expressed  the  mean  velocity  as  flow  rate 
per  unit  area.  For  both  samples,  the  military  specifications  call  for  50 
to  90  cfm/ft2  at  a  differential  pressure  of  .5  inch  of  water,  and  450  to 
650  cfm/ft2  at  a  differential  pressure  of  20  inches  of  water.  Sample  B 
extrapolated  to  88  cfm/ft2  and  654  cfm/ft2,  respectively,  and  sample  C 
extrapolated  to  89  cfm/ft2  and  487  cfm/ft2,  respectively. 

Figure  4  shows  the  pressure  coefficient  Cp,  defined  by 
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Figure  3.    Differential  Pressure  vs.  Dynamic  Head 
for  Parachute  Cloth  Samples 


10 


U">. 
CD 


o 
1\ 


LO. 

U3 


3 

to 


LJ 

O 

L_  lo 
L_  uV 
LJ 
O 
O 

LJ  o_ 

to 
en 

LJ 

rv-    LT> 


\ 


V 


o 

LEGEND 
1.1   02.  NYLON 

A 
+ 
X 

4.75  0Z.  NYLON 

7.0  0Z.  NYLON 

'""3:0""0Z"."KEVLRR 

0.00         0.01  0.02  0.03  0.04 

DYNAMIC  HERD 


0.05 


0.06 


0.07 


Figure  4.   Pressure  Coefficient  vs.  Dynamic  Head 
for  Parachute  Cloth  Samples 
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for  all  the  samples  as  a  function  of  dynamic  head.  Clearly,  the  pressure 
coefficient  remains  nearly  independent  of  velocity  for  a  given  canopy. 
As  noted  earlier,  the  flow  on  the  downstream  face  of  the  sample 
(which  ranged  from  0  to  14.4  ft/sec)  had  no  effect  on  the  results 
shown  in  Figures  3  and  4.  Table  2  shows  the  minimum,  maximum, 
and  average  values  of  the  pressure  coefficient  for  the  different 
samples. 

TABLE  2 

COMPARISON  OF  THE  PRESSURE  COEFFICIENTS  FOR 
THE  DIFFERENT  CANOPY  MATERIAL  SAMPLES 

Sample      Minimum  Cp       Maximum  Cp  Average  Cp 

A  20.3  27.3   (32.4)  22.9   (24.5) 

B  33.0  43.0   (52.0)  36.7   (39.2) 

C  39.0  48.9   (54.8)  42.8   (44.8) 

D 61.5 73.3   (80.4)  66.3   (69.1) 

Notes: 

( 1 )  The  values  in  parentheses  in  the  maximum  column 
are  for  the  lowest  velocity  data  points. 

(2)  The  average  values  do  not  include  the  pressure  coeffi- 
cients for  the  lowest  velocity  data  points  because  they 
were  deemed  not  representative  of  the  velocities  that 
parachutes  descend.  The  values  in  parentheses  in  the 
average  column  are  averages  if  these  data  points  had 
been  included. 
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HI.    NUMERICAL   EXPERIMENTS 

A.      THEORETICAL  FOUNDATIONS  OF  THE  MODEL 

The  flow  about  a  cambered  plate  is  obtained  by  successive  trans- 
formations from  a  circle  plane.  In  fact,  it  can  be  shown  that  the  flow 
in  the  £  -plane  (circle  plane)  may  be  transformed  to  the  flow  in  the 
z-plane  (camber  plane)  through  the  use  of 


z  =  £+  m  - 


£  +  m 


2m2-  1 
m 


in  which  m  is  a  parameter  related  to  the  total  angle  of  the  circular  arc. 
For  example,  m  equal  to  .5,  .707,  or  .866  corresponds  to  camber 
angles  of  120,  180,  or  240  degrees,  respectively.  Figure  5  shows  the 
transformation  from  a  circle  to  a  120-degree  camber. 


Figure  5.    Circle  and  Physical  Planes 
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The  complex  function  describing  the  flow  is  given  by 


+  ^Ln 


in  which  Tkp  and  £kp  represent  respectively  the  strength  and  location 
of  the  k-th  p- vortex  (counter-clockwise  rotating),  and  Tkq  and  £kq 
represent  the  strength  and  location  of  the  k-th  q-vortex  (clockwise 
rotating).  As  described  in  Munz,  discrete  vortices  are  introduced  at 
the  tips  of  the  camber  at  regular  time  intervals,  and  evolution  of  the 
vortex  wake  is  calculated  through  the  use  of  the  velocities  induced  at 
each  vortex  location.  Furthermore,  the  forces  acting  on  the  camber 
are  calculated  through  the  use  of  the  generalized  Blasius  theorem,  or 
through  the  integration  of  the  differential  pressure.  The  numerical 
model  requires  the  input  of  a  velocity  history,  and  in  the  present 
investigation,  the  velocity  profile  that  was  employed  was  that  used  by 
Munz  for  the  simple  purpose  of  delineating  the  effect  of  porosity 
relative  to  the  non-porous  case  calculated  by  Munz. 

The  existing  analysis  had  to  be  considerably  modified  in  order  to 
incorporate  the  effect  of  porosity.  The  inclusion  of  the  effect  of 
porosity  was  already  shown  to  be  difficult  as  far  as  the  previous 
experimental  investigations  are  concerned.  It  turns  out  that  its  math- 
ematical analysis  is  just  as  difficult  as  its  experimental  understanding. 
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B.     INOUES  METHOD 

The  only  work  which  dealt  with  a  steady  uniform  flow  past  a  sim- 
ple porous  flat  plate  was  carried  out  by  Inoue.  To  emphasize,  the  flow 
was  not  subjected  to  any  deceleration  and  the  only  unsteadiness  intro- 
duced into  the  flow  came  from  the  evolution  of  the  vortex  shedding 
following  the  impulsive  start  of  the  flow  at  a  constant  velocity.  Inoue 
made  no  attempt  to  relate  the  differential  pressure  across  the  plate  to 
the  flow  rate.  According  to  his  method,  the  complex  function  for  a 
plate  normal  to  the  flow  is  given  by 


f  = 


a(C-t)      (l-a)(C  +  t) 


i2  2 

h  r      t-c       vrL     C-CL 

XKi     *     ^K      •  Z,— —In ir- 
7-ln- — r+1      2ji      c-I 

in  which  a  is  an  artificial  porosity  coefficient  related  to  the  geometric 
porosity  and  aUmf  is  the  part  of  the  flow  which  passes  through  the 
plate  uniformly,  regardless  of  the  differential  pressure  distribution 
prevailing  along  the  plate.  Inoue  did  not  calculate  forces  or  pressure 
distributions,  and  confined  himself  to  the  comparison  of  the  calculated 
flow  kinematics  with  sample  flow  pictures,  presumably  obtained  under 
conditions  similar  to  those  encountered  in  the  numerical  calculations. 
Inoue's  analysis,  though  not  applicable  to  the  analysis  of  the  actual 
flow  under  consideration,  was  tried  anyway,  to  gain  some  experience 
regarding  the  overall  performance  of  the  code  developed  in  the  pres- 
ent investigation.  Figures  6  through  10  show,  at  representative  times, 
the  evolution  of  flow  for  a  =  0.25  and  for  the  velocity  profile  given  by 
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Figure  6.    Wake  Vortex  Positions  for  a  =  0.25  at  T  =  8.85 
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Figure  7.   Wake  Vortex  Positions  for  a  =  0.25  at  T  =  11.29 
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Figure  8.   Wake  Vortex  Positions  for  a  =  0.25  at  T  =  12.92 
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Figure  9.    Wake  Vortex  Positions  for  a  =  0.25  at  T  =  14.54 
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Figure  10.    Wake  Vortex  Positions  for  a  =  0.25  at  T  =  16.98 
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U/Uq  =  1  for  T  =  U0t/c  <  9.72 
and 

U/Uq  =  0.97T  -  0.05T2  -  3.70 
and 

A  =  0.97  -0.10T 
in  the  interval  9.72  <  T  <  11.48  and 

U/Uq  =  -0.3423  +  0.0072T2  +  3.82 
and 

A  =  -0.3423  +  0.01445T 

in  the  interval  11.48  <  T  <  17.95. 

For  T  greater  than  18,  the  velocity  and  acceleration  are  zero. 

Figures  11  through  15  show  the  flow  at  exactly  the  same  times 
about  a  nonporous  camber  (a  =  0).  Comparison  of  these  figures  shows 
that,  as  anticipated, the  effect  of  the  porosity  is  to  push  the  wake  fur- 
ther downstream  and  thereby  minimize  the  consequences  of  the 
return  of  the  wake  to  the  body.  In  fact,  Figures  16  through  20  show 
the  comparison  of  the  differential  pressures  at  the  corresponding 
times  for  a  =  0.0  and  0.25.  Apparently,  the  differential  pressure  for  a 
=  .25  becomes  less  negative  than  the  one  for  a  =  0.0,  even  though,  as 
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Figure  11.    Wake  Vortex  Positions  for 
Non-Porous  Camber  at  T  =  8.85 
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Figure  12.   Wake  Vortex  Positions  for 
Non-Porous  Camber  at  T  =  11.29 
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Figure  13.   Wake  Vortex  Positions  for 
Non-Porous  Camber  at  T  ■  12.92 
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Figure  14.   Wake  Vortex  Positions  for 
Non-Porous  Camber  at  T  =  14.54 


9- 


25 


CO 

en 

LD 

• 

to 

CD 

• 

*— ' 

oo 

1 

1 

O 

D_ 

,.  .    .    / 


:*. 


E- 


V- 


i 


^  i 


•n 
I 


t- 


ZI 


Figure  15.   Wake  Vortex  Positions  for 
Non-Porous  Camber  at  T  =  16.98 


26 


-60        -50         -40  -30         -20 


-10  0  10 

THETR 


40  50  60 


Nonporous  Camber 


1 1 1 1 1 ^-t 1 1 1 1 r 

-60        -50         -40  -30  -20  -10  0  10  20  30  40  50  60 

THETR 

Porous  Camber,  a  =  0.25 
Figure  16.    Differential  Pressure  Comparison  at  T  =  8.85 
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Porous  Camber,  a  =  0.25 
Figure  17.    Differential  Pressure  Comparison  at  T  =  11.29 
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Porous  Camber,  a  =  0.25 
Figure  18.    Differential  Pressure  Comparison  at  T  =  12.92 
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Porous  Camber,  a  =  0.25 
Figure  19.    Differential  Pressure  Comparison  at  T  =  14.52 
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Figure  20.    Differential  Pressure  Comparison  at  T  =  16.98 
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noted  earlier,  Inoue's  method  does  not  take  into  account  the 
relationship  between  the  flow  and  the  differential  pressure,  or  the 
variability  of  a  along  the  camber  and  with  time.  In  other  words,  a 
must  vary  with  time  and  with  the  angle  along  the  camber.  Thus,  the 
comparison  of  the  a  =  0.0  and  a  =  0.25  cases  will  not  be  made  in 
greater  detail  than  Inoue's  method  deserves. 

C.   PRESENT  METHOD  OF  POROSITY  MODELING 

With  the  realization  of  the  fact  that  the  flow  through  the  camber 
must  be  related  to  the  prevailing  differential  pressure  and  that  the 
direction  of  the  flow  in  the  wake  has  no  measurable  effect  on  the  dif- 
ferential pressure  versus  dynamic  head,  a  new  method  had  to  be 
devised  which  will  properly  account  for  the  effect  of  porosity. 

The  present  method  takes  advantage  of  the  simplicity  of  Inoue's 
method  and  accounts  for  the  local  dependence  of  the  flow  on 
differential  pressure.  For  this  purpose,  the  value  of  a  is  chosen  so  as 
to  make  the  Vn  at  0  =  0.0  proportional  to  the  V  ^P  as  given  by  Figure 
4.  The  variations  in  the  flow  across  the  camber  beyond  and  above  that 
given  by  Inoue's  method  are  accounted  for  through  the  use  of  discrete 
vortices  placed  along  the  camber.  A  number  of  sample  calculations 
have  shown  that  the  use  of  118  body  vortices  (one  every  degree)  on  a 
120-degree  camber  will  be  more  than  sufficient  to  account  for  the 
additional  cross-flow  effects  smoothly.  The  first  vortices  were  placed 
one  degree  from  the  tips,  and  no  vortex  was  placed  at  6  =  0.  The  mid- 
points  between   the  vortices  were   used   as   the   control  points   to 


32 


calculate  the  velocity  normal  to  the  camber.    At  each  control  point,  the 
differential  pressure  was  calculated,  and  then  Vn  was  found  from 


v  =J-/^E 
N    cpv    p 

The  velocity  induced  by  the  uniform  flow  based  on  a  was  subtracted, 
and  then  the  strengths  of  the  body  vortices  were  determined  by 
equating  the  remaining  velocity  to  that  induced  by  the  body  vortices  at 
the  control  points.  Once  the  strengths  of  the  body  vortices  were 
determined,  the  velocity  induced  at  each  vortex  in  the  wake  was  cal- 
culated, the  velocity  of  the  ambient  flow  was  updated,  and  all  of  the 
wake  vortices  were  convected  for  a  time  interval  At  using  a  2nd  order 
modified  Eulerian  scheme.  Then  the  procedure  was  repeated,  up  to 
any  specific  time  (in  this  analysis  T  =  24).  The  wake  vortices  were 
never  merged  together,  no  matter  how  close  they  came,  and  were  only 
removed  from  the  simulation  if  they  contacted  the  camber.  To  pre- 
vent abnormally  high  induced  velocities  from  occurring  due  to  the 
proximity  of  one  vortex  to  another,  whether  body  or  wake  vortices, 
Rosenhead's  (1930)  smoothing  scheme  was  used,  where  the  complex 
velocity  at  a  point  was  multiplied  by 

x2 


x2+52 


where  x  is  the  distance  from  the  point  to  the  vortex  inducing  the 
velocity  and  5  is  the  attenuation  factor.  In  this  analysis,  8  was  chosen 
to  be  0.05. 
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The  method  described  above  is  quite  general  and  could  be  used  to 
investigate  the  effect  of  the  porosity  of  any  canopy.  Because  it  is 
deemed  sufficient  to  understand  the  physics  of  the  phenomenon,  only 
two  canopy  characteristics  were  used  in  the  numerical  experiments. 
The  first  case  corresponded  to  the  characteristics  of  canopy  sample  A 
(Cp  =  22.9,  see  Table  2),  and  the  second  to  a  canopy  which  is  much 
more  porous  than  any  of  the  samples  tested,  with  an  assumed  Cp  = 
4.45,  which  corresponded  to  an  initial  value  of  a  =  0.25  at  9  =  0  and 
T  =  .5. 

Figures  21  through  25  show  the  evolution  of  the  wake  for  Cp  = 
22.9  (sample  A).  Comparison  of  these  figures  with  Figures  11  through 
15  shows  that  the  effect  of  porosity,  however  noticeable,  is  not  large 
enough  to  show  a  significant  downstream  displacement  of  the  wake 
vortices.  The  importance  of  this  conclusion  lies  in  the  fact  that  the 
most  porous  parachute  canopy  used  by  the  industry  is  not  sufficiently 
porous  to  minimize  the  consequences  of  the  wake  return,  i.e.,  the 
collapse  of  the  parachute  as  evidenced  by  field  experiments.  Thus,  it 
has  been  possible  through  a  relatively  idealized  analysis  to  show  the 
inadequacies  of  the  materials  used  in  the  construction  of  large 
parachutes. 

The  next  step  in  the  analysis  was  to  show  the  beneficial  effects  of 
porosity  when  it  can  be  made  sufficiently  large.  For  this  purpose,  a 
value  of  Cp  =  4.45  was  chosen,  for  the  reasons  noted  earlier.  Figures 
26  through  30  show  the  evolution  of  the  wake  and  Figures  3 1  through 
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Figure  21.   Wake  Vortex  Positions  for  Cp  =  22.9  at  T  =  8.85 
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Figure  22.    Wake  Vortex  Positions  for  Cp  =  22.9  at  T  =  11.29 
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Figure  23.   Wake  Vortex  Positions  for  Cp  =  22.9  at  T  =  12.92 
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Figure  24.    Wake  Vortex  Positions  for  Cp  ■  22.9  at  T  =  14.54 
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Figure  25.    Wake  Vortex  Positions  for  Cp  =  22.9  at  T  =  16.98 
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Figure  26.    Wake  Vortex  Positions  for  Cp  =  4.45  at  T  =  8.85 
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Figure  27.   Wake  Vortex  Positions  for  Cp  =  4.45  at  T  =  11.29 
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Figure  28.    Wake  Vortex  Positions  for  Cp  =  4.45  at  T  =  12.92 
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Figure  29.   Wake  Vortex  Positions  for  Cp  =  4.45  at  T  =  14.54 
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Figure  30.    Wake  Vortex  Positions  for  Cp  =  4.45  at  T  =  16.98 
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35  the  differential  pressures.  It  is  evident  from  a  comparison  of  Fig- 
ures 26  through  30  with  either  Figures  11  through  15  (nonporous 
case)  or  Figure  21  through  25  (Cp  =  22.9  case)  that  the  wake  is 
significantly  displaced  downstream  and  that,  as  seen  in  Figures  31 
through  35,  the  differential  pressure  is  far  less  negative  than  those 
shown  in  Figures  16  through  20  (a  =  .25  and  a  =  0.0).  Figure  36 
shows  histories  of  the  differential  pressure  computations  at  0  =  0 
degrees  for  the  nonporous  case,  the  present  model  with  Cp  =  22.9, 
Inoue's  method  with  a  =  0.25,  and  the  present  model  with  Cp  =  4.45. 
It  is  easily  seen  that  the  present  method  with  Cp  =  22.9  has  differen- 
tial pressures  only  slightly  less  in  magnitude  than  the  nonporous  case, 
while  the  two  other,  more  porous  cases  have  differential  pressures 
significantly  smaller  in  magnitude  and  delay  the  occurrence  of  nega- 
tive differential  pressure  at  the  center  of  the  body.  It  is  also  seen  that 
while  the  two  large  porosity  cases  during  steady  state  and  the  early 
deceleration  periods  are  very  similar,  during  later  stages  of  the  decel- 
eration the  present  method  results  in  lower  differential  pressures  as  it 
accounts  for  the  amount  of  flow  through  the  body  due  to  the  local  dif- 
ferential pressures. 
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Figure  31.    Differential  Pressure  for  Cp  =  4.45  at  T  ■  8.85 
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Figure  32.    Differential  Pressure  for  Cp  =  4.45  at  T  =  11.29 


Figure  33.    Differential  Pressure  for  Cp  =  4.45  at  T  =  12.92 
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Figure  34.    Differential  Pressure  for  Cp  =  4.45  at  T  =  14.54 


Figure  35.    Differential  Pressure  for  Cp  =  4.45  at  T  =  16.98 
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Figure  36.    Differential  Pressure  History  Comparison 
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IV.    CONCLUSIONS 

The     results     presented     herein     warranted     the     following 
conclusions. 

1.  The  effect  of  porosity  is  a  function  of  the  differential  pressure 
across  the  porous  material,  and  is  dictated  only  by  the  physical 
characteristics  of  the  material. 

2.  For  a  given  canopy,  the  normalized  differential  pressure  (the 
pressure  coefficient)  is  independent  of  the  ambient  velocity,  up 
to  a  velocity  of  about  600  miles  per  hour. 

3.  The  effect  of  porosity  cannot  be  adequately  accounted  for  by  let- 
ting through  the  camber  a  constant  uniform  velocity,  as  done  by 
Inoue.  The  effect  of  the  variation  of  differential  pressure  with 
time  and  along  the  camber  must  be  included. 

4.  A  novel  method  has  been  devised  which  combines  the  uniform 
porous  flow  concept  with  the  variability  of  the  flow  rate  with  dif- 
ferential pressure  through  the  use  of  discrete  vortices  embedded 
in  the  camber. 

5.  The  results  have  shown  that,  in  general,  the  effect  of  porosity  is 
to  displace  the  wake  vortices  further  downstream  and  to  retard 
their  return  to  the  camber  during  the  period  of  deceleration. 

6.  Even  the  most  porous  canopy  material  used  in  the  construction  of 
large  parachutes  is  not  sufficiently  porous  to  prevent  the  collapse 
of  the  parachute,  as  evidenced  by  the  field  experiments  and  veri- 
fied by  the  present  numerical  model. 

7.  The  results  of  calculations  with  a  much  larger  porosity  have 
shown  that  it  will  be  possible  to  considerably  delay  or  entirely 
eliminate  the  collapse  phenomenon.  In  doing  so,  it  may  be 
advantageous  to  optimize  the  radial  distribution  of  porosity  along 
the  camber  in  order  to  achieve  as  large  a  drag  force  as  possible. 
One  must,  however,  bear  in  mind  the  fact  that  one  also  needs  to 
account  for  the  flexibility  of  the  canopy.  This  forms  the  basis  of 
the  continuing  efforts  toward  the  understanding  of  the  behavior 
of  aerodynamic  deceleration  devices. 
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